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I. INTRODUCTION 

Many of the fascinating properties of correlated elec¬ 
tron systems derive from collective behavior. In systems 
with a magnetic impurity embedded in a host of con¬ 
duction electrons, the Kondo effecfPl represents a prime 
example. It is caused by the antiferromagnetic exchange 
interaction between the impurity spin and the local spin 
of conduction electrons. Below a characteristic tempera¬ 
ture Tk, the Kondo temperature, the magnetic moment 
of the impurity is collectively screened by a cloud of con¬ 
duction electrons with spatial extension determined by 
Ck vf/Tk^ where vf is the Fermi velocity.!^ Kondo 
screening is a phenomenon that eludes a perturbative ap¬ 
proach, as is expressed by the non-analytic dependence 
of Tk on the strength J of the local exchange coupling: 
IuTk ^ —IjJ. For weak J, the Kondo scale is exponen¬ 
tially small. This translates into a mesoscopically large 
Kondo cloud. 

Hence, for a magnetic impurity embedded in a small 
nanostructure®^ thn concept of Kondo screening must 
be modified.l^^n^®®^ The universal aspects of Kondo 
physics are lost when Tk is of the same order of mag¬ 
nitude as the finite-size gap A at the Fermi energy of 
the conduction-electron system or, if viewed from a real- 
space perspective, when the Kondo cloud exceeds the size 
of the system. A Kondo scale in the regime of some tens 
of Kelvins roughly corresponds to a linear extension in 
the nanometer range. 

If there are several magnetic impuritiesthe sit¬ 
uation gets complicated by the competition of Kondo 
screening with non-local magnetic correlations. The lat¬ 
ter are typically induced by an indirect exchange cou¬ 
pling which is mediated by the conduction electrons. 
The effective Ruderman-Kittel-Kasuya-Yosida (RKKY) 
inter act iorP^M^ can be derived perturbatively and scales 


as Jrkky ^ In the weak-J regime this domi¬ 

nates over the exponentially small Kondo scale, while the 
Kondo effect sets in for sufficiently strong J > 

This competition between Kondo screening and RKKY 
interaction is also qualitatively modified by quantum 
confinement For coupling strengths weaker than Ja, 
where Ja is defined by the condition A = Tk{Ja)^ 
the conventional Kondo effect is replaced by a “finite- 
size Kondo effect” where an impurity spin is screened 
by forming a singlet with a single conduction electron 
occupying the fully delocalized one-particle state at the 
Fermi energy. The energy to break this singlet scales lin¬ 
early (rather than exponentially) with J and therefore 
the finite-size Kondo effect dominates over the RKKY 
interaction for J -^0. Hence, with decreasing J there is 
a reentrance of the Kondo effect 

In this paper we focus on the conditions under which 
different magnetic moments (impurity spins) in a nanos¬ 
tructure can be individually screened. To investigate pos¬ 
sible screening mechanisms on the linear-in-J scale, we 
consider a generic multi-impurity Kondo model involv¬ 
ing a conduction-electron system with a finite number of 
sites. In the weak-coupling limit J < Ja, the physics 
of this model is accessible by standard perturbation the¬ 
ory. In particular, we study the case where the electronic 
structure of the non-interacting conduction-electron sys¬ 
tem exhibits more than one single-particle state at the 
Fermi energy sf due to symmetries of the confinement 
geometry. There are several questions to be answered in 
this context, e.g.: Which states at Sf provide a “chan¬ 
nel” for the screening of an impurity spin? How many 
impurity spins can be screened for J ^ 0? In which 
way does the geometry and the electronic structure of 
the nanosystem affect the screening? 

Gapless Kondo models (A = 0) with higher spin quan¬ 
tum numbers or multiple screening channels^^^^^ have 
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attracted considerable interest in the past since their 
low-energy p ropert ies cannot be captured by local Fermi- 
liquid theoryP^l^ The present paper addresses the sim¬ 
pler but analogous question: in a nanosystem with sev¬ 
eral magnetic impurities and degenerate orbitals at the 
Fermi energy, when and how does an overall spin-singlet 
ground state arise? Are there analogs of underscreening 
and overscreening? 

A highly interesting question concerns the crossover 
from the physics of the “finite-size multi-channel Kondo 
effect” to the physics of the gapless system, realized in the 
thermodynamic limit. This crossover should be visible 
when extending the nanostructure in size, corresponding 
to the limit A ^ 0 and Ja ^ 0 , or when increasing J 
and thereby approaching the regime J > Ja where finite- 
size effects no longer matter. With increasing coupling 
strength, the individual screening on the linear-in-J scale 
is expected to be replaced by non-local magnetic corre¬ 
lations due to the operation of the RKKY coupling at 
order until, on further increasing the coupling toward 
the strong-J limit, the standard Kondo effect is again 
recovered. 

To elucidate our ideas with concrete examples we 
consider the two-impurity Kondo model on a one¬ 
dimensional lattice in a ring geometry. This model not 
only shows individual screening of the two impurity spins 
for J ^0, depending on details of the geometry and elec¬ 
tronic structure, but also exhibits various crossovers as 
a function of J and of the system size. These are stud¬ 
ied numerically by employ ing th e density-matrix renor¬ 
malization group (DMRG]P3^ as we ll as the numeri¬ 
cal renormalization group (NRG) We use a matrix- 
product implementation and a folded-chain geometry to 
account for the periodic boundary conditions within the 
DMRG. Application of the NRG implementation to an 
impurity (or multi-impurity) system with a discrete bath 
requires a specific adaption of the theory and the numer¬ 
ical algorithm to perform the Wilson-chain mappings. 

The paper is organized as follows: After briefly in¬ 
troducing the model in Sec. [HI the general perturbative 


theory is explained in Sec. |III| Results obtained with the 
DMRG and the NRG for the two-impurity Kondo model 
on a finite ring are discussed in Secs. and |V| respec¬ 
tively. A summary and the conclusions are given in Sec. 

ED 


II. MULTI-IMPURITY KONDO MODEL 

Using standard notations, the multi-impurity Kondo 
model on a lattice of finite size is given by the following 
Hamiltonian: 

R 

H = Hq -f- H\ = ^ ^ (7 T J ^ ^ ^ir ‘ • (1) 

ii'a r=l 

The first part Hq describes a system of N non-interacting 
conduction electrons moving on a lattice of arbitrary di¬ 
mension and geometry consisting of L sites. There is a 


spin-degenerate orbital |i,cr) at each site i, and a =t, i 
is the spin projection, cj^ and Qcr are the correspond¬ 
ing creators and annihilators. Since we wish to address 
the physics of magnetic nanostructures, we take L to 
be finite. Diagonalization of the spin-independent tight- 
binding hopping matrix ta ^, 

=£khk'^gg' , ( 2 ) 

ii' 

yields the set of one-particle energies £k and one-particle 
eigenstates |/i:,^,cr), where the index g = l,...,G(/c) ac¬ 
counts for the possible degeneracy of Sk arising, e.g., due 
to spatial symmetries of the nanostructure. The uni¬ 
tary matrix U has elements Ui-kg = and is 

spin-independent. Three-dimensional systems with lin¬ 
ear extension in the range of ^ 10 nm correspond to a 
level spacing Ack = £k-\-i — of the order of some tens 
of Kelvins. 

The second part Hi describes R magnetic impurities, 
modeled as spins with spin quantum number S = 1 / 2 . 
Each impurity is coupled by a local antiferromagnetic 
exchange J > 0 to the conduction-electron system. The 
r-th spin Sr interacts with the local conduction-electron 
spin density Si^ = | at site Here, 

a denotes the vector of Pauli matrices. 

The physics of this model in the thermodynamic limit 
with 1 / ^ oc is hi ghly complex and characterized by 
a subtle competitiorP^^^ between Kondo screening on 
a scale given by the KondcP temperature Tk and the 
indirect RKKY^f^ magnetic coupling Jrkky- While 
the RKKY coupling is a perturbative concept and ob¬ 
tained with second-order-in-J perturbation theory, the 
Kondo effect is non-perturbative with InTx ^ — 1/J. Be¬ 
cause the Kondo temperature is exponentially small, the 
RKKY coupling Jrkky ^ J^ dominates in the weak-J 
regime. 

As was pointed out in Ref. |24l this picture may change 
qualitatively for systems of finite size L when Tk ^ A, 
where A is the finite-size gap at the Fermi energy. In 
this case, the logarithmic Kondo correlations are cut by 
the finite system size,^ and the universal Kondo effect 
is replaced by a “finite-size Kondo effect”. The latter is 
characterized by a linear-in-J energy scale, whose pre¬ 
cise value depends on details of the geometry and elec¬ 
tronic structure, and which dominates over the RKKY 
interaction as J ^ 0. Since the finite-size gap A regular¬ 
izes the problem, the “finite-size Kondo effect” should be 
well accessible to standard perturbation theory for cou¬ 
pling strengths J < Ja, where Ja is given by the condi¬ 
tion A = Tk(Ja)- Since InT^ ^ “1/J foi* Ihe infinite 
system,[Ilwe have J < — 1 /lnA. 


III. EFFECTIVE LOW-ENERGY THEORY 

Let |FS, 7 ), with 7 = l,...,r, be the (F-fold degener¬ 
ate) Y-electron ground state of Hq (the Fermi sea) which 
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is obtained by occupying all one-particle levels below the 
Fermi energy, Sk < 

One can distinguish between two situations: In the 
“off-resonance” case, the ground state is non-degenerate, 
r = 1, and the orbitals \kY^g^ <j) with ^ = 1,..., G(/cf) at 
= ^/cf fully occupied. Note, that this requires an 
even number of electrons in the nanostructure. In this 
off-resonance case, the interaction Hi can be treated as 
a perturbation to recover the effective RKKY model. 



rr' 


where JRKKY,rr' is given by times the non-local 
(r 7 ^ r') static spin susceptibility of the non-interacting 
conduction-electron system. In particular, there is no 
linear-in-J Kondo effect in this case. 

In the following, we will concentrate on the “on- 
resonance” case where F > 1, i.e., where the one-particle 
orbitals with energy sy are incompletely occupied. Note 
in particular that this is necessarily the case when the 
nanostructure hosts an odd number of electrons N. Em¬ 
ploying degenerate perturbation theor}^!! up to linear or¬ 
der in J, the low-energy physics is captured by an effec¬ 
tive model 


= P^HiPn , (4) 

where Pq = 1^ u projector onto the 

space of ground states of the unperturbed (J = 0) Hamil¬ 
tonian. To compute PoSi^Po, we first consider the uni¬ 
tary transformation 

^\gcr ~ , (5) 

i 

which gives 

^ir — 2 Pkg]ir^\ga^cr(^' g' a'Uir-,k' g' • ( 6 ) 

kk' ,gg' ,aa' 

Since the Pauli matrices are traceless, terms with k^k' ^ 
k^ do not contribute to PoS^^Pq: 

Po^i^Po = - PkFg]ir^l^Fgcr^^^'^^Fg'(T'Pir]kFg' Po • 

gg',aa' 

„ ( 7 ) 

The effective low-energy model Eq. Q can be written as 
a spin-only model. To this end we introduce another set 
of site-dependent unitary transformations: 

^kFOiA'^r) = '^^^ga{'i^r)<^kFgo- ' (^) 

9 

In general, the unitary G(/cf) x G(/cf) matrices V{ir) 
are different for each site v (and do not commute); 
the creator is therefore also dependent on p. 

V(ir) are the transformations that diagonalize the dyadic 


products Ugg>{ir) = Ui^-kFg' for each with r = 

1,...,P: 

^agi'^r)Ul^a.i^Uir]kFg'^g'a' (v) — Xa(ir)^aa' • (9) 

99 ' 

Trivially, this construction implies that there is a sin¬ 
gle non-zero eigenvalue, Xa=i(ir) = 'E^g \Ui^-,kF 9 \‘^• The 
corresponding eigenvector has the elements Vga=i{ir) = 

^Lg-,ir/\j^g IF,;fepgP- Using this, and collecting the 
results, we find 

PoSi^Po = Po ^ \Ui^-kFg\‘^^F{ir) , (10) 
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where 


^f(v) — 2 ^kpaai'^r)o'aa' 


C/CF 


/ (i^) 


a=l 


( 11 ) 


is the conduction-electron spin on the (spin-degenerate) 
Eermi orbital |E,P,cr) = ^(iF)|vac.). Explicitly, 

we have 

|F> U, cr) = / == E E ^r) • 

y Eg \k'ir-,kFg\ g i 

( 12 ) 

The goal of the transformation ® is to cast the effec¬ 
tive Hamiltonian in a spin-only form. With the definition 
of the effective coupling. 


Teff(v) = J ^ \Uir-Mg?‘ 5 (13) 

9 

Eq. Q follows as 


R 


HqH — Pq ^ ^ Teff (^r) ^F (“^r) ’ 


(14) 


r=l 


Eormally, the effective Hamiltonian has the structure of 
a central-spin model with the subspace of conduction- 
electron states at the Eermi level as the “central” degrees 
of freedom. However, the low-energy physics at weak J 
decisively depends on different orbitals |E,i^,cr) that de¬ 
fine the conduction-electron spins Sf(v) lo which the 
impurity spins Sr are coupled. These orbitals are degen¬ 
erate one-particle eigenstates of Hq with energy which 
are delocalized over the entire lattice. Since |E,P,crWs 
in general different for each impurity spin Sr, Eq. (14) 
may represent an unconventional central-spin model. 

Depending on the details of the geometry and the elec¬ 
tronic structure, the values of Jeff (v) and also the defini¬ 
tions of Sf(v) fo Eq- (14) can be very different. Several 
distinct situations can arise, which we discuss now in 
turn. 

Eirstly, the coupling constant Jeff(ir) foi* a particular 
impurity spin Sr can vanish if |Pv;/cf^P = 0 fo^* all g 
(i.e., each conduction-electron wave function Ui-kFg l^as 
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a node at the site v). In this case, there is no finite-size 
Kondo screening of the impurity at site ir on the linear- 
in-J scale, and second-order perturbation theory includ¬ 
ing RKKY couplings must be employed to understand 
the lowest-energy physics. By contrast, an impurity at 
another site ir' with finite Jeff(V') can still be screened 
by forming a delocalized spin-singlet state with a con¬ 
duction electron at This situation was discussed in 
Ref. [24] for the case with G(/cf) = 1 but T = 2 (a sin¬ 
gle level at the Fermi energy, with an odd total number 
of electrons), in the context of a one-dimensional Kondo 
model at half-filling with up to three impurities. 

Secondly, we consider the case G(/cf) = 1 and odd 
Y, as above, but for a Fermi wave function = 

Ui^-kFg=i which is non-vanishing at all impurity sites v. 
An illustration is given by Fig. [^(a) and (b). Depending 
on details of the geometry, the resulting coupling Jeff(ir) 
is possibly site-dependent. Since G(/cf) = 1, the matrices 
V{ir) reduce to 1 x 1 “matrices” trivially {V{ir) = 1 = 
1), and all impurity spins Sr couple to the spin of the 
same orbital |F,cr) = |F, v,cr). In this case, the effective 
model reduces to a conventional central-spin model, 

R 

i^eff = PqSf • ^ Jeff (ir)‘^r , (15) 

r=l 

ii i2 

Si S2 


A • 



FIG. 1: (Color online) Schematic illustration of the ef¬ 

fective low-energy theory: (a) Multi-impurity Kondo model 
with hopping t and antiferromagnetic exchange coupling J. 
Example with two impurity spins Si and S 2 coupled to 
the conduction-electron spins at sites ii and Z 2 of a one¬ 
dimensional chain. (b) One-particle eigenenergies Sk for 
J = 0. A single electron occupies the non-degenerate one- 
particle state at the Fermi energy sf- A (A') is the finite-size 
gap. For 0 < J t the effective low-energy model i/eff 
(dashed line) is given by a central-spin model with effective 
couplings Jeff (A) between Sr and the spin sf of the delocal¬ 
ized one-particle state at sf- Only a single screening channel 
is available in this case. 


with a single central spin sf which is independent of 

This situation can be viewed as an antiferromagnetic 
Heisenberg model on a bipartite “lattice”, with the cen¬ 
tral spin Sf on an “A sublattice” site and the impurity 
spins Sr on the “B sublattice” sites. Hence, the Lieb- 
Mattis theorenP^ applies: The ground state has total 
spin 5'tot = {P ~ l)/2, i-e., there is exactly one screen¬ 
ing channel. For R > 2, this results in “underscreening” 
of the impurity degrees of freedom (whereby the over¬ 
all ground state remains degenerate). The intuitive pic¬ 
ture is that all impurity spins are antiferromagnetically 
aligned to the central spin, (sf • < 0, which implies 

ferromagnetic impurity spin correlations {Sr • Sr') > 0. 

This underscreening phenomenon can be seen most 
clearly in the special case of a translationally-invariant 
D-dimensional finite lattice comprising L < oo sites (for 
example, a ring geometry in D = 1). The periodic bound¬ 
ary conditions in this case imply that the modulus square 
of the Fermi wave function Ui^-kFg independent of the 
site index ir. Hence 


Deff = To JeffSF • ^ Sr 


(16) 


The model can then be straightforwardly solved by us¬ 
ing the standard rules for adding angular momenta: The 
ground state of iJeff in Eq. (16) with ground-state energy 


Eeff,o = — Jeff(T + l)/4 is obtained for maximum total 
impurity-spin quantum number 5'imp = ^/2 and for total 
spin quantum number 5'tot = Simp “ consistent with 
the Lieb-Mattis theorem, and is (2S'tot + l)-fold degen¬ 
erate. The ^-components of Sf and of the total impu¬ 
rity spin iS'imp = '^r Pr in the ground state with max¬ 
imal magnetic quantum number Mtot = S'tot are given 
by ~ T ^/-R +1 and {SiYnp,z) ~ 

respectively. 

Finally, and most interesting, we discuss the case 
where each impurity spin Sr couples to a different 
conduction-electron spin, corresponding to the Fermi or¬ 
bital |F, i^,cr). The inner product of the Fermi orbitals 
is given by 


(F, V,cr'|F, V,cr) = 


Pir''^^FgPkFg-,i 


\j^g \l^9 


(17) 

Suppose that all orbitals |F,i^,cr) are mutually orthogo¬ 
nal and that the number of impurity spins R is equal to 
the degeneracy of the Fermi energy G{kY). In this case, 
the effective model (14) describes a Heisenberg model 


which decomposes into a set of R = G(/cf) decoupled an¬ 
tiferromagnetic subsystems, each consisting of an impu¬ 
rity spin Sr and the corresponding conduction-electron 
spin s{ir). The ground state is then a total spin sin¬ 
glet — there are as many screening channels as impurity 
spins, and one has “perfect” screening. The degeneracy 
of the one-particle eigenenergies at the Fermi level, and 
mutual orthogonality of the states |F, D, a), typically re- 
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suits from spatial symmetries of the system. A concrete 
example of this is discussed in Sec. [IV| 

The case R > G(/cf) corresponds to underscreening, 
since the number of available orthogonal Fermi orbitals is 
insufficient to form independent singlets with all impurity 
spins. Consequently, some of the impurity spins do not 
couple to the conduction-electron systems (to first order 
in J), leaving a spin-degenerate ground state. 

On the other hand, if R < G(/cf), there are more 
screening channels than necessary. However, each impu¬ 
rity spin Sr is coupled to at most one conduction-electron 
orbital |F,v,cr), precluding the frustrated “overscreen¬ 
ing” of any given impurity. In this case, some of the 
Fermi electrons simply remain decoupled from the impu¬ 
rities, giving rise to a trivial ground-state degeneracy. 

If the inner product 0 is finite, the situation is more 
complicated because the creators and annihilators defin¬ 


ing the conduction-electron spins in Eq. (14) do not refer 
to orthogonal orbitals. Generally this leads to an un¬ 
conventional central-“spin” model, where the Fermi elec¬ 
trons can interact with more than one impurity. This 
case requires a different numerical approach, and will be 
discussed in a separate paper.l^ 


IV. CROSSOVERS BETWEEN STRONG 
COUPLING AND KONDO OR RKKY REGIMES 

Having explored the expected physics of finite-sized 
magnetic nanostructures on general grounds, we turn 
now to a specific example with R = 2 impurity spins 
coupled to a conduction-electron system with doubly- 
degenerate orbitals at the Fermi energy (G(/cf) = 2). 

Going beyond the analysis of the previous section, here 
we also consider the various crossovers expected to arise 
upon increasing the strength of the exchange coupling, 
J. These crossovers are manifest in the evolution of zero- 
temperature static spin-spin correlation functions, which 
we calculate numerically. In particular, when the two im¬ 
purities are locked together into a spin-singlet, we expect 
{Si ' S 2 ) = — for the inter-impurity spin correlation. 
This situation pertains when the impurities are coupled 
by an effective non-local antiferromagnetic RKKY inter¬ 
action. On the other hand, if the Kondo effect dominates 
(and provided both impurity spins can be screened), we 
expect {Sr • Stot) = for the correlation between each 
impurity spin r = 1,2 and the total conduction-electron 
spin Stot = Si. The inter-impurity correlation should 
vanish in this case, (iS'i • S 2 ) =0, signaling that the im¬ 
purities are decoupled from each other. 

The question of which process dominates is a subtle 
one, since the RKKY coupling cannot be treated as an 
independent parameter but rather depends on the Kondo 
coupling J. Furthermore, besides the RKKY coupling 
-^RKKY and the Kondo temperature Tk, the finite-size 
gap A represents a third energy scale. The competition 
between these is the focus of our present study. 

As a concrete example, we now consider a one- 



FIG. 2: (Color online) Left: Two-impurity Kondo model 
on a finite one-dimensional lattice with periodic boundary 
conditions. Off-resonance case with L — 4m -h 2 (illustrated 
for m = 1). Right: Schematic picture of the one-particle 
eigenenergies Ek for J = 0. Eigenenergies (except for the 
lowest and highest) are two-fold degenerate, G(/cf) = 2. ep 
denotes the highest occupied one-particle energy. At half¬ 
filling N — L the one-particle states at are fully occupied. 


dimensional lattice with a finite, but possibly large num¬ 
ber of sites L, to which two spin-^ impurities are cou¬ 
pled. This two-impurity Kondo model can be solved for 
arbitrary coupling strength J by numerical means, as dis¬ 
cussed further below. 

Specifically, we assume a non-zero hopping tij = —t 
between nearest neighbors i and j of a one-dimensional 
lattice, and set t = 1 to fix the energy scale. Two impu¬ 
rity spins iS'i and S 2 are coupled to nearest neighbor sites 
ii and 22 - We realize two-fold degenerate Eermi orbitals, 
G(/cf) = 2, by imposing periodic boundary conditions 
(i.e. a ring geometry). The model is studied at zero tem¬ 
perature and half-filling {N = L electrons). 

The hopping matrix is diagonalized [c.f. Eq. @1 by 
discrete Eourier transformation, 

Ui,kg = > ( 18 ) 


with the “-h” sign for ^ = 1 and “—” sign for ^ = 2, wave 
vector k = —(tt/u) -\-nAk (where n = 0,1,..., L — 1), and 
position vector Ri = i• a for site i = 0,1,— 1. Here, a 
is the lattice constant and Ak = 27r/(al/). Note that the 
effective couplings (13) do not vanish since \Ui.kg\‘^ 7^ 0- 
The two-fold degenerate one-particle energies are given 
by 


ek = —2tcos{ka) . (19) 

This results in a finite-size gap at the Eermi energy of 
A = 2t sin(27r/I/) ex 1/L for large L. 

Eor moderate system sizes up to L = 100 and interme¬ 
diate coupling strengths, numerical results for the static 
T = 0 spin-spin correlators (iS'i • S 2 ) and {Sr • Stot) can 
be accurately obtained usiM the density-matrix renor¬ 
malization group (DMRG).^ We implement here a stan¬ 
dard scheme based on matrix-product states and matrix- 
product operators; see Ref.[42]for a brief discussion of the 
algorithm, and Ref. |36| for a general overview of matrix- 
product-state techniques. 

Usually, DMRG is formulated for and applied to one¬ 
dimensional systems with open boundaries. To employ 
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FIG. 3: (Color online) Evolution of the spin-spin correlation 
functions {Si • ^ 2 ) (top) and {Si • Stot) (bottom) with J, 
as obtained from DMRG for various system sizes L. Off- 
resonance cases for L = 4m+ 2 (lines), on-resonance cases for 
L = 4m with integer m (dots). See Fig. for a schematic of 
the setup. 


the standard algorithm for the present case with peri¬ 
odic boundaries, we first consider the open chain consist¬ 
ing of L sites and connected by nearest-neighbor hopping 
terms and for i = 1, ...,1/ — 1 as usual. The 

chain is folded in half and connected at the open end 
with the missing hopping terms ti^L and tL,i to gener¬ 
ate a new half-length chain with new “sites” consisting 
of pairs (i,I/ — i + 1) of original sites. We let the stan¬ 
dard DMRG algorithm operate on the new “sites”. The 
approach avoids long-range hopping terms in the Hamil¬ 
tonian (which are unfavorable for the scaling of DMRG), 
at the cost of an enlarged local Hilbert space. We address 
larger system sizes and weaker J in Sec. |V| using an alter¬ 
native approach based on the numerical renormalization 


We discuss first off-resonance situations, realized at 
half-filling for L = 4m + 2 with integer m — see Fig. 

for an illustration. Here the ground state for J = 0 is 
the non-degenerate Fermi sea, and all one-particle states 
at sy are completely filled. Recall from Sec. |m] that no 
linear-in-J finite-size Kondo effect can arise in the off- 
resonance case, because all Fermi orbitals are completely 


occupied. Instead, the effective low-energy physics in the 
limit J ^ 0 is described by the RKKY model, Eq. (§. 
The sign of the RKKY exchange coupling depends on 
the distance d = \ii —i 2 \ between the impurity spins; for 
the present case of a half-filled one-dimensional system, 
-^RKKY cx: (—1)^+^. In particular, for nearest neighbors 
{d = 1) the coupling is antiferromagnetic^ ^/rkky > 0- 

Fig. shows DMRG results for the inter-impurity cor¬ 
relations {Si • S 2 ) and for the correlation between one 
of the impurity spins and the total conduction-electron 
spin {Si • Stot)- Note that reflection symmetry implies 
{Si • Stot) = {S 2 • Slot). Results have been obtained for 
different system sizes L = 10, 50, 98. Fig. [^demonstrates 
that the impurities are indeed effectively decoupled from 
the conduction-electron system for J ^ 0, as is indicated 
by the vanishing correlation function {Si • Stot) ^ 0. 
The antiferromagnetic RKKY exchange, coupling impu¬ 
rity spins Si and S 2 , is dominant at weak coupling. This 
results in the formation of a non-local RKKY singlet, 
with {Si • S 2 ) as J ^ 0. 

With increasing J one observes a smooth crossover 
from this RKKY regime to a strong-coupling regime 
where the two impurity spins are separately screened by 
the conduction electrons. This is reflected by a vanishing 
inter-impurity correlation { 81 - 82 ) ^ 0 as J ^ oc, and 
by Kondo correlations (^S^-Stot) ^ For large J, the 
Kondo effect simply reduces to the formation of entirely 
local singlets, i.e., we find {Sr ' Si^) — 2/4 for the local 
Kondo correlation function (not shown in the figure). 

The system-size dependence of the crossover from 
RKKY to Kondo physics is quite regular in the off- 
resonance case, and the results for L = 98 are rather 
representative of the physics in the thermodynamic limit 
L ^ 00 . This is demonstrated explicitly and discussed 
further in Sec. m The off-resonance behavior here can be 
understood from the conventional Doniach scenarioit 
is associated with a characteristic coupling strength Jd 
separating RKKY and Kondo regimes. 

However, the on-resonance situation (realized here for 
L = 4m, see Fig. [^ is more complex, since there one 
expects a second crossover on the scale of Ja to a regime 
dominated by the linear-in-J finite-size Kondo effect. We 
explore this physics now in the context of full DMRG 
results with L = 12,48,100 presented in Fig. 

First, we comment briefly on the strong-coupling limit 
J ^ Td, Ta- Here the physics is essentially the same as 
that arising in the off-resonance case, both being domi¬ 
nated by strong local spin correlations {Sr • Details 
of the electronic structure close to the Fermi level then 
become irrelevant. This can be seen directly by compar¬ 
ing the off- and on-resonance situations (lines vs. dots) 
in Fig. For example, the same basic behavior is found 
for all J > 1 with L = 98 and L = 100 (green lines and 
dots); and agreement becomes essentially quantitative for 
J > 3 or so. 

By contrast, on- and off-resonance cases are qualita¬ 
tively very different at weak coupling. Here, the low en¬ 
ergy details of the conduction-electron system become 
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FIG. 4: (Color online) Same as Fig. but for an on- 

resonance case with L = 4m (illustrated here for m = 2). 
At half-filling N — the one-particle states at sf are par¬ 
tially occupied. 


important. In the on-resonance case, the one-particle 
states at the Fermi energy are half-filled, as sketched 
in Fig. The resulting degeneracy of the many-body 
ground state at J = 0 gives rise to a linear-in-J Kondo 
scale, which dominates as J ^ 0. The effective low- 
energy model is then given by Eq. (14). 

Indeed, our DMRG results are indicative of a vanish¬ 
ing inter-impurity correlation (iS'i • iS' 2 ) ^ 0 as J ^ 0, 
implying that the RKKY interaction is inoperative at 
weak coupling. Strong Kondo correlations do however 
appear to develop, with {Sr • Stot) ^ as J —^ 0. 
This is seen clearly from the numerics with L = 12 (and 
results for L = 48 are highly suggestive). However, we 
were unable to stabilize well-converged DMRG results for 
larger systems and smaller J, due to the need to resolve 
small energy scales A ^ l/L, and the development of 
highly extended entangled ground states. Nevertheless, 
the results are consistent with our expectation that the 
ground state for L = 4m is an overall spin-singlet, with 
each impurity being separately screened. 

In contrast to the strong-coupling (large J) limit, 
where the screening can be visualized as two spatially 
separated and almost perfectly localized Kondo singlets, 
for small J < Ja the “Kondo clouds” extend over the 
entire system (i.e. they have extent L). This implies 
the existence of two different screening channels. Simi¬ 
larly, in the thermodynamic limit L ^ oc, it was pointed 



FIG. 5: (Golor online) Schematic picture of the effective 
low-energy model of the system depicted in Fig. The site- 
dependent unitary transformation V(ir) generates orthogonal 
Fermi orbitals. 


out in Ref. [22] that perfect and essentially independent 
Kondo screening can take place for nearby impurities, 
even though their Kondo clouds are extended and sub¬ 
stantially overlapping in real-space. The difference, how¬ 
ever, is that as 1/ ^ 00 the conduction-electron system 
becomes gapless A ^ 0, and the Kondo cloud for a single 
impurity then has extent ^ '^f/Tk (with up the Fermi 
velocity), which can be understood from a renormaliza¬ 
tion group perspective.!^ 

The DMRG results are fully consistent with the impli¬ 
cations of the effective low-energy Hamiltonian (14). For 


J < Ja the two impurities interact with a central region 
given by the one-particle states at and the two Fermi 
electrons. To determine whether the impurities couple to 
orthogonal states or not, we evaluate Eq. 0 using the 
wave functions given by Eq. (H^: 


( 20 ) 


(F,ii,(T'|F,i 2 ,cr) = cos [fcF(-Rii - -Ris)] 


Note the different dependence on the distance as com¬ 
pared to the RKKY coupling, which is Jrkky oc 
cos[2/cp(i7q — Ri^)]- At half-filling, the Eermi wave vec¬ 
tors are k-p = 4z7rl2a. With \Ri — R 2 I = a • d, we 
have (E, cr|E, 22 , cr) = cos(7rd/2). In the case of near¬ 
est neighboring impurities (d = 1) the matrix element 
vanishes, implying that the impurity spins couple to or¬ 
thogonal channels, and are screened separately (see Eig. 

To summarize these findings: Eor very large J, each 
impurity is Kondo-screened by conduction-electron states 
which are highly localized in real-space. As J is de¬ 
creased, the Kondo screening becomes more delocalized. 
As J is reduced further, below Jd, there is a crossover 
to an inter-impurity RKKY singlet state. In the off- 
resonance case, the RKKY interaction continues to dom¬ 
inate as J ^ 0. But in the on-resonance case, a second 
crossover arises on the scale of Ja to a state charac¬ 
terized by the finite-size Kondo effect, with orthogonal 
Kondo clouds spreading over the entire lattice. Provided 
Ja do, these regimes will be distinct, and re-entrant 
Kondo physics should be seen. 

Since A ^ l/L in one dimension and thus Ja ^ 
l/lnL, one generally expects that two successive 
crossovers will be cleanly observed upon decreasing J for 
larger system sizes. This is of course physically sensi¬ 
ble, since in the L ^ 00 limit, the second (finite-size 
Kondo) crossover is pushed to J ^ 0, and becomes un¬ 
observable. By contrast, for small L (see e.g. DMRG 
results for L = 12 and 48 in Eig. [^, Ja and Jd are not 
well-separated. Then, the RKKY regime is never fully 
realized, and (iS'i • S 2 ) > — for all J. 

Einally, we touch on the physics arising when the im¬ 
purity spins are separated by an even distance, d. Here, 
the modulus of the overlap (20) is unity, implying that 


both impurity spins couple to the spin in the same Eermi 
orbital. The ground state in this case is a tensor product 
of the (underscreened) doublet ground state of the cor¬ 
responding central-spin model, and the Kramers doublet 
of the Eermi sea (with the remaining unpaired electron 
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V. APPROACHING THE THERMODYNAMIC 
LIMIT 

For substantially larger systems, one expects that 
Ja ^ Jt>‘ In this case, a well-defined RKKY regime 
should develop, manifest by inter-impurity spin correla¬ 
tions approaching {Si • S' 2 ) ^ —3/4. Furthermore, the 
distinction between on- and off-resonance cases must be¬ 
come irrelevant in the thermodynamic limit L ^ 00 for 
coupling strengths J ^ Ja ^ l/lnL ^ 0. Then, the 
conduction-electron system becomes gapless, and it is 
unimportant that states precisely at the Fermi energy 
are either fully or partially occupied. 

As discussed in Sec. El larger systems and weaker cou¬ 
plings cannot in practice be investigated with DMRG. To 
understand the evolution on increasing the system size 
to approach the thermodynamic limit, we now employ 
a different but complementary approach, based on the 
numerical renormalization groupP^ (NRG). 

Traditionally, NRG is used for single impurities cou¬ 
pled to gapless conduction-electron systems in the ther¬ 
modynamic limit (for a review, see Ref. [38|) . In that case, 
NRG is known to provide accurate and efficient a ccess 
to a wide range of thermodynamic!^ and dynamicP^Eil 
quantities, at essentially any temperature or energy 
scale. Although generalization of NRG to systems 
with several conduction-electron channels requ ires sig¬ 
nificantly greater computational resour ces,l^^^^ detailed 
results for two impurities separated in real-space on a 
three-dimensional lattice with L 00 were recently ob¬ 
tained in Ref. |2^ Note that NRG can be used for impu¬ 
rity models of any geometry and dimensionality, unlike 
DMRG which is usually confined to one spatial dimen¬ 
sion. 

To study finite systems with NRG, we generalize the 
standard formulation. In particular, two features of the 
current problem are non-standard and must be treated 
differently. Firstly, the conduction-electron density of 
states for the finite system consists of a set of discrete 
poles located at Sk^ rather than a continuum. Secondly, 
in the on-resonance case of particular interest, there is 
a pole located at the Fermi energy, which controls the 
low-energy physics. At the heart of NRG lies a logarith¬ 
mic discretization of the conduction-electron spectrum; 
in the present context, the spectrum is already discrete, 
and a pole at the Fermi energy cannot be resolved on 
the logarithmic scale. Our NRG implementation, which 
overcomes these issues, is discussed in the Appendix. 

In the following we discuss our NRG results, obtained 
for the effective two-impurity, two-channel setup using 
a logarithmic discretization parameter A = 2 and re¬ 
taining Nk = 15000 states at each step of the itera¬ 
tive process. We exploit total charge and spin-projection 
quantum numbers (but so-called z-trick averaging was 
not used). In all cases, we have verified that the results 


are fully converged with respect to indicating that 
the discretized model is solved exactly by NRG. 

Fig. shows NRG results (lines) for the full J- 
dependence of the inter-impurity spin correlation {Si • 
iS' 2 ), obtained for different system sizes from L = 10 up 
to 1/ = 10^. For small L, the previously discussed DMRG 
results are included for comparison (circle points). Off- 
resonance cases {L = 4m + 2) are shown in black; on- 
resonance {L = 4m) in blue. 

For very large system sizes, both on-resonance {L = 
10^) and off-resonance {L = 10^ + 2) cases agree per¬ 
fectly with results in the thermodynamic limit L ^ 00 
(red cross points), where NRG should give numerically 
exact results.l^ In the other extreme, for very small sys¬ 
tem sizes 1/ = 10 and L = 50 (off-resonance) or L = 12 
and L = 48 (on-resonance), NRG agrees very well with 
the available DMRG data. This is a non-trivial result, 
since NRG involves a re-discretization of the conduction- 
electron spectrum, which itself only contains I//2 poles 
(see Appendix). However, for the on-resonance case 
L = 100, we start to see a more significant deviation 
between the depth of the minimum in (iS'i • S 2 ) calcu¬ 
lated using DMRG and NRG. This discrepancy can be 
traced to the discretization in NRG. Results obtained us¬ 
ing different A indicate that this discrepancy diminishes 
as A ^ 1, where the bare model is recovered. We also 
found for a given A > 1 that the DMRG results are recov¬ 
ered using NRG with a slightly larger L, i.e., the effect 
of discretization in NRG can be also be thought of as 
a slight renormalization of the system size L. The gen¬ 
eral trend exhibited as L increases is, however, consistent 
with expectation from Secs. |III| and [IV| as now discussed. 


We find that the results for the off-resonance case L = 
4m+ 2 converge uniformly and rapidly with increasing L. 
In all such cases, (iS'i • S 2 ) ^ —3/4 as J —> 0, indicating 
the dominant RKKY mechanism at weak-coupling. 

In the on-resonance cases with finite but large L = 
1000 and L = 10000, the spin correlations {Si • S 2 ) 
displayed in Fig. clearly reflect a dominant RKKY 
coupling mechanism in the intermediate-J range. In 
particular, the minimum in (iS'i • S 2 ) approaches —3/4, 
symptomatic of inter-impurity singlet formation. In this 
regime, the corresponding value of {Sr • Stot) = ~| ~ 
(iSi •iS' 2 ) ~ 0 indicates vanishing Kondo correlations. The 
crossover from the strong-coupling Kondo to this RKKY 
regime is nearly independent of system size and is located 
at Jd ^1 — 2. 

On the other hand, the crossover from the RKKY 
regime to the finite-size two-channel Kondo screened 
ground state takes place on the scale of Ja ^ 1/ In L. For 
J ^ Ja, the spin-spin correlations coincide with the con¬ 
tinuum result {L —> 00 ); only in the range 0 < J < Ja 
does one see marked deviations. However, this regime 
progressively shrinks to zero as Ja ^ 1 /lnL ^ 0, indi¬ 
cating that the weak-coupling limit J ^ 0 and the ther¬ 
modynamic limit L ^ 00 do not commute. For system 
sizes L = 100 and smaller, Ja and Jd are of compara¬ 
ble size; then only incipient RKKY effects are observed. 
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FIG. 6: (Color online) Evolution of {S\ • S^) with J, as obtained from NRG for various system sizes L (on-resonance cases: 
blue lines, off-resonance cases: black lines). Crosses: continuum result in the limit L ^ oo obtained by standard NRG. DMRG 
data from Fig. are shown for comparison (circles). 


The clean crossover to a distinct RKKY regime is only 
observed for larger L. 


The non-trivial weak-coupling physics is further ex¬ 
amined in Fig. The lower panels show the charac¬ 
teristic energy scale which is necessary to break the 
singlet ground state for both off-resonance (left panel) 
and on-resonance (right panel) cases. This energy scale 
also controls the renormalization-group flow to the sta¬ 
ble fixed point describing the ground state within the 
NRG scheme. Clearly, one would expect that Tg oc 
in the off-resonance case corresponding to the effective 
RKKY model (§ and Tg oc J in the on-resonance case 
corresponding to the central-spin model (IT3. To check 
this, we have calculated with NRG the full (zero temper¬ 
ature) dynamical correlation function {{Si; S 2 ))uj^ i-e., 
the Fourier transform of iO{t){Si{t) ■ >S'2(0)), following 
Ref. im As cj ^ 0, we find {{Si; S 2 ))uj -> 0 - the van¬ 
ishing of the correlator indicates spin-singlet formation. 
At finite uj ^ Tg^ however, S2))uj shows a peak re¬ 
sulting from spin-flip fluctuations. Hence, Tg can be in¬ 
terpreted as the energy scale for singlet formation. From 
the lower panels of Fig. [^we can read off the following 


asymptotic behavior: 

Tg (off resonance) , 

Tg J (on resonance) , (21) 


which nicely agrees with our expectations, and confirms 
the perturbative analysis. 

In the upper panels of Fig. the asymptotic scaling 
behavior of (S^i •iS' 2 ) is analyzed for the weak-J limit. For 
the off- and the on-resonance case we find 


(^1 • ^2) 





(off resonance) , 


(^1 • ^2) 


j^o 



(on resonance) , 


( 22 ) 


respectively. Here, J = J^o + a/hi{hL) has a weak loga¬ 
rithmic dependence on system size, and Lq ~ 300 (a and 
b are constants). Within the perturbative regime defined 
by Eqs. and ( p!4| ) , the spin correlations cannot exhibit 
a J dependence (there is only a single ener gy s cale). The 
corrections oc and oc J given by Eq. (22) therefore 
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FIG. 7: (Color online) Upper panels: J dependence of {Si ■ S 2 ) for small J for various system sizes L. Lower panels: J 
dependence of the low-energy scale Tg, corresponding to singlet formation. Distinct asymptotic scaling is found for off-resonance 
cases (left panels) and on-resonance cases (right panels). Results obtained by NRG. 


reflect the J-dependence of terms in higher-order pertur¬ 
bation theory. 


VI. CONCLUSIONS 

The physics of multiple magnetic impurities embed¬ 
ded in a host metal is characterized by the competition 
between Kondo screening and RKKY interactions. In 
quantum confined nanostructures, the picture changes 
qualitatively and the physics is far richer. In particular, 
there are parity effects, and a new regime can emerge at 
weak coupling dominated by finite-size effects. 

For a single impurity, the finite-size effect 
can be viewed in two complementary ways. If the 
conduction-electron spectrum has a gap A around the 
Fermi energy that is larger than the Kondo temperature 
Tk , then spin-flip scattering and the Kondo effect is sup¬ 
pressed. If there is no conduction-electron state at the 
Fermi energy ep (the so-called off-resonance case), the 
impurity remains unscreened down to T = 0. But in 
the on-resonance case, a state precisely at ep can bind 
with the impurity and screen its spin through singlet for¬ 
mation. This is the finite-size Kondo effect, and scales 
linearly in the impurity-host coupling J. Since the finite- 
size gap cuts off the Kondo correlations on the scale of A, 
the physics on the lowest energy scales is perturbatively 
accessible. 

Alternatively, one can think of the problem in real- 


space: if the Kondo screening cloud of extension can¬ 
not fit inside a nanostructure of finite size L, then the 
Kondo effect is replaced by its finite-size variant, in which 
impurity screening involves a delocalized conduction- 
electron state spanning the entire system. 

In finite systems with several impurities,l^^^^ the 
RKKY scale must also be considered. In the off- 
resonance case, the impurities become coupled by the 
RKKY interaction on the lowest energy scales. For two 
impurities, they can collectively either form a singlet or 
a triplet state, depending on whether the RKKY interac¬ 
tion is antiferromagnetic or ferromagnetic. Which is real¬ 
ized in practice depends on details of the nanostructure 
geometry and electronic structure, and is encompassed 
by standard RKKY theory. 

But, since the RKKY interaction scales as J^, the 
finite-size Kondo effect dominates at weak coupling in 
the on-resonance case. This physical picture implies 
that for a particular system, on decreasing the coupling 
J, one can realize two successive crossovers: first from 
a standard Kondo strong-coupling state to the RKKY- 
dominated regime, and then on further decreasing J, to 
the finite-size Kondo regime. At weak coupling J ^ 0, 
the finite-size Kondo effect always wins out in the on- 
resonance case. But if several impurities are present, can 
they all be screened by finite-size Kondo effects? If not, 
does underscreening lead to degenerate ground states? Is 
frustrated overscreening a possibility? 

In this paper we answer these questions by putting 
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the general concepts on firm footing. We focused mainly 
on the on-resonance case, deriving first an effective gen¬ 
eralized central-spin model at weak coupling. All of the 
impurities (despite being located at different sites in real- 
space) couple to the same small number of conduction- 
electron states which lie precisely dXep- Several scenarios 
can arise, depending on the nanostructure geometry and 
electronic structure. We show that it is entirely possible 
for each impurity to be coupled to its own Fermi orbital, 
allowing exact finite-size-Kondo screening of all impuri¬ 
ties, and an overall spin-singlet ground state. In other 
situations, there are insufficient degrees of freedom dX sp 
to screen all of the impurities. The remaining unscreened 
impurities are then coupled together via their mutual 
RKKY interaction. Depending on details, this could ei¬ 
ther lead to a degenerate ground state or a singlet ground 
state, with both finite-size-Kondo and RKKY mecha¬ 
nisms acting simultaneously on different impurities. We 
show that overscreening can never occur, since each im¬ 
purity is coupled to at most one conduction-electron state 

dX Ep. 

The implications of the effective low-energy effective 
model were substantiated by full DMRG and NRG calcu¬ 
lations for a system comprising two impurities on a finite 
one-dimensional ring of conduction-electron sites. This 
model, although simplified, captures a number of features 
relevant to magnetic nanostructures. In particular, de¬ 
pending on the system size, one can realize either on- or 
off-resonance cases. Furthermore, the coupling between 
the impurities and the Fermi orbitals in the on-resonance 
case can realize either exactly screened or underscreened 
scenarios, depending on impurity separation. 

We also investigated numerically the full crossovers 
between finite-size-Kondo, RKKY, and regular Kondo 
strong-coupling physics on increasing J for a system of 
given size L, going beyond the perturbative analysis. For 
smaller systems, where DMRG can be used, the RKKY 
regime is never fully realized, since there is competition 
between RKKY and finite-size-Kondo effects. But for 
larger systems, accessible with NRG, the two crossovers 
are cleanly observed, with a distinct intermediate RKKY 
regime. The final crossover to finite-size-Kondo is pushed 
to weaker and weaker coupling as L increases. Approach¬ 
ing the thermodynamic limit L ^ oc, this crossover is 
pushed to J ^ 0, as intuitively expected. 

It would be interesting to explore the physics of more 
realistic finite three-dimensional systems with several 
impurities. Analysis of the low-energy effective the¬ 
ory in each case should elucidate the screening mecha¬ 
nisms. Full numerical calculati ons coul d be implemented, 
using a block-Lanczos methocP^^^^^^ to map the real- 
space system onto a chain (or ladder) form amenable to 
treatment with either DMRG or NRG. Indeed, ffatband 
models,E^M^ which sustain a macroscopically large num¬ 
ber of states at the Fermi energy, could also be investi¬ 
gated. 

Finally, an interesting open question regards the fate 
of the impurities when the conduction-electron states at 


ep are not perfectly degenerate, as would arise due to 
symmetry-breaking perturbations (e.g. real-space distor¬ 
tions or coupling anisotropies). Then one would expect 
different screening mechanisms to be available at differ¬ 
ent energy scales. In this context, also the broadening of 
the energy levels due to a residual coupling of the nanos¬ 
tructure to the environment must be considered.l^^K^ An¬ 
other related question concerns particle-hole symmetry 
breaking, and the relation to hard-gapped problems,!^ 
where quantum-phase transitions between Kondo and 
unscreened states might also be accessible. 
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Appendix A: NRG calculations 


In this App endix , we describe our implementation of 
Wilson’s NRG,I^^^^ generalized to treat a system com¬ 
prising two spin- 1/2 impurities coupled to neighboring 
sites of a fermionic tight-binding ring of L = 2 m sites. 
The physical setup is illustrated in Fig. [^a). 

The conduction-electron Hamiltonian, Hq in Eq. 0- 
can be written in the form. 


m —1 


Ha = 


E (< 


n=l cr,a=± 

+‘E( 


ancr^a(n+l)cr 


+ H 


.c.) 


(Al) 


^+lcr^-lcr + ^+mcr^-mcr 


H.c. 


The impurity part of the Hamiltonian is Hi = J(S+i • 
s+i-fS_i-s_i). We now transform to an even/odd orbital 
basis for the conduction electrons, defined in terms of the 
canonical fermions, C(e/o)na = :^(c+na ± c_na), viz: 


m —1 

^o=iE E G 


n—1 (j^oc—elo 


.t r 

ancT Q/(n-|-l)cr 


+ H.C.) 






C.l, - c 




olcT^olcT ' ^ema^ema ^oma^oma 


(A2) 


Eor J = 0, the even and odd channels are strictly decou¬ 
pled. Eor finite J, each impurity couples to both even and 
odd channels — see Eig.|^b). This two-impurity problem 
is thus manifestly and irreducibly two-channel in nature. 
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FIG. 8: (Color online) Mappings f or th e two-impurity model 
in NRG. (a) Real-space basis Eq. (Al), with two impurities 
coupled to neighboring sites of a tight-binding ring comprising 
a hnite even number of fermionic sites L = 2m. (b) Even/odd 
basis Eq. (A2), with impurities coupled to both even and odd 
conduction-electron ch anne ls, each comprising L/2 sites, (c) 
Discretized model Eq. (A6), with impurities coupled to even 
and odd Wilson chains, each comprising N < L/2 sites. 


To use NRG, we must map each channel a = e/o into 
the form of a Wilson chain.123 

NRG is usually used for gapless Fermi systems, and 
involves a logarithmic discretization of the (continuous) 
conduction-electron density. In the present context of 
the finite ring, the J = 0 conduction-electron density of 
states pa{^) for each channel a = e/o is already discrete, 

L 

Pa{oj) = '^aa{p)S [w - e{p)] . (A3) 

p=l 

However, the poles located at s{p) are not of course 
distributed on a logarithmic energy grid, as required 
for NRG. Employing the diagonal representation Eqs. 
(§, lU), ([T^ (together with the even/odd basis trans¬ 
formation), we find e{p) = 2tcos[27r£] and a^/oip) = 
^(1 ±cos[27r£]). 

To perform the Wilson chain mappings, we now re¬ 
discretize (or re-bin) each spectrum pa{^) on a loga¬ 
rithmic frequency grid. This is done by dividing pa{^) 
into intervals of exponentially-decreasing width, defined 
by the discretization points = ±2tA“’^ (with n = 
0,1,2,3... and A > 1). All poles in a given interval are 
replaced by a single pole of the same total weight to yield, 

= E ±)5 [w - Un. ±)] , (A4) 

n,± 

where, 

L 

ba{n, ±) = E aa{p)0[±e{p) T x^_^_^]e[±x^ =F £{p)] , 

p=l 

Un,±) = E T ^t+im±xi T e{p)] ■ 

(A5) 

For a given real-space system with L orbitals, the number 
and distribution of poles in p^^®‘^(c(;) is thus controlled by 


off-resonant: on-resonant: 



FIG. 9: (Golor online) Wilson chain coefficients x 
and x vs Wilson index n using K — 2 for various 

system sizes L. 


the discretization parameter A. In general there are (far) 
fewer poles in p^^®^(cc;) than in the true pa{^)- 

For e ach ch annel a, the Wilson chain is now defined 
uniquel}E^^ as the one-dimensional tight-binding chain 
which has the same local density of states at one end as 
the discretized spectrum p^^®^(co’). 


N 

^o"^ = E[E 


a,cr n=0 
N-1 


(A6) 


E^“ (/fn./„(„+l). + H.C.) 


n=0 


where N = N{L^ A), = e^(I/, A) and = t^{L^ A) de¬ 

pend on the original system size L and the discretization 
parameter A. The number of Wilson chain orbitals N 
in each channel is finite and typically (far) smaller than 
the original number of real-space orbitals L. In prac¬ 
tice, the Wilson chain mapping is achieved by Lanczos 
tridiagonalization,!^^^^ usin g as input the pole weights 
and positions from Eq. (A5). The impurity subsystem is 
then coupled to the end of the Wilson chains, as depicted 
in Fig.j^c). 
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1. Fermi-level pole 


As discussed in Sec. the underlying physics of the 
model is expected to be very different, depending on 
whether or not the conduction-electron system has a pole 
at the Fermi level. With system size L = 2m as above, 
no Fermi level pole exists in the off-resonance case for 
odd m, while a pole of weight 2/L lies precisely at the 
Fermi level in the on-resonance case for even m. 

When there is no Fermi level pole, we can safely use 
the logarithmic discretization scheme embodied by Eq. 
(A5) to generate the Wilson chain formulation, Fig.j^c). 

However, in the on-resonance case, the important ef¬ 
fect of the Fermi level pole cannot be captured by a dis¬ 
cretization scheme defined on a logarithmic grid. To this 
end, we broaden the Fermi level pole us ing a Lorentzian 
of width 6 . Specifically, we replace Eq. (A3) by 


However, the results converge rapidly with decreasing 
5. Finally, we use S = 10“^^. Although physically this 
pole broadening should have negligible effect (being far 
smaller than any physical energy scale of the problem, 
including T^), its inclusion does affect the Wilson chain 
coefficients, allowing the physics of the on-resonance case 
to be studied with NRG. 

Fig. shows the resulting Wilson chain coefficients for 
the same systems as in Fig.|^ Note that for both on- and 
off-resonance cases the Wilson chain is of finite length 
N (meaning that = 0). Their non-trivial evolution 
encodes the real-space physics of the tight-binding ring. 


2. Iterative diagonalization 




L 

+ ^ X 


S/t 


p=i 




(A7) 


where the primed summation indicates that the pole at 
e(p) =0 is excluded. The discretized pole positions and 
weights from Eq. (A5) are therefore modified to include 


the extra density in each interval. 


ba{n, ±) =ba{n^ i) ^ — arctan 
27r 




^n+iy 


LS 

ia{n, ±)ba{n, ±) =Ca{n, ±)ba{n, ±) ± — In 


(A8) 


Note that any finite broadening 6 results in an infinite 
number of discretized poles (albeit with exponentially de¬ 
creasing weight). The Wilson chain mapping can now be 
performed as before, and the results examined as a func¬ 
tion of 6 . 

We find that inclusion of the Fermi level pole has a sig¬ 
nificant effect on the resulting Wilson chain coefficients. 


The key property of the Wilson chain hoppings 
is that they decay exponentially down the chain — 
Fig.ll This is obviously not a property of the phys¬ 
ical real-space system, but arises due to the logarith¬ 
mic (re-)dis cretiz ation. As with Wilson’s original NRG 
formulationj ^^ l ^^ l this justifies an iterative process of di¬ 
agonalization and truncation. 

Starting from the impurity subsystem, one builds up 
the chains by successively coupling on additional Wilson 
orbitals. At a given step n < N, corresponding to a sys¬ 
tem comprising the impurities and n Wilson orbitals of 
each channel, the Hamiltonian is diagonalized. However, 
only the lowest-energy Nk states are kept for construc¬ 
tion of the Hamiltonian at step n + 1. Because decay 
exponentially down the chain, the discarded high-energy 
states at iteration n do not cross over into the low-energy 
manifold at a later iteration n' > n. The physics of the 
system is therefore examined on successively lower energy 
scales as more Wilson orbitals are added, embodying the 
RG structure of the problem. 

The spin-spin correlators presented in Fig. were 
obtained from the full thermal density mat rixp ^ con¬ 
structed in the complete Anders-Schiller basis.ESl 
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